Doped Mott Insulators are Insulators: Hole localization in the Cuprates 
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We demonstrate that a Mott insulator lightly doped with holes is still an insulator at low tem- 
perature even without disorder. Hole localization obtains because the chemical potential lies in a 
pseudogap which has a vanishing density of states at zero temperature. The energy scale for the 
pseudogap is set by the nearest-neighbour singlet-triplet splitting. As this energy scale vanishes if 
transitions, virtual or otherwise, to the upper Hubbard band are not permitted, the fundamental 
length scale in the pseudogap regime is the average distance between doubly occupied sites. Con- 
sequently, the pseudogap is tied to the non-commutativity of the two limits U — ► oo (U the on-site 
Coulomb repulsion) and L — > oo (the system size). 



Hole doping a Mott insulator shifts the chemical 
potential from the middle of the charge gap generated 
by the energy cost (17) for double occupancy to the top 
of the lower Hubbard band. Nominally, the density of 
states at the top of the lower Hubbard band is non- 
zero. Consequently, doped Mott insulators are expected 
to be conductors. However, doped Mott insulators such 
as the high temperature cuprate superconductors are well 
known to possess a pseudogap 00 d at the Fermi energy 
below some characteristic temperature, T*, that persists 
well into the superconducting dome. While a dip in the 
density of states is not sufficient to destroy the simple pic- 
ture that a metallic state obtains upon light hole doping, 
certainly a vanishing density of the states at the Fermi 
level would be. The question arises: Does the density of 
states vanish at the chemical potential in the limit T — > 
in the underdoped cuprates or in lightly doped Mott in- 
sulators in general? The analysis presented here on the 
Hubbard model suggests the answer to this question is 
yes and hence lightly doped Mott insulators are, in fact, 
still insulators. 

Experimental probes that shed light, either directly or 
indirectly, on the ultimate fate of the density of states 
at the chemical potential in doped Mott systems are of 
three types: transport, tunneling and angle-resolved pho- 
toemission (ARPES). Early transport measurements on 
La 2 _ a; Sr a; Cu04 + j / revealed^ @ that in the lightly-hole 
doped regime, the in-plane resistivity obeys the 3-d vari- 
able hopping form, 



p(T) oc e^^° 
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with a = 1/40 or diverges logarithmically as lnTo/T0. 
In both cases, if localization is due to disorder (and 
hence extrinsic to Mott physics), an externally applied 
magnetic field should couple to the orbital motion and 
yield a negative magnetoresistance. While the mag- 
netoresistance is negative, it is independent of the di- 
rection of the fieldgg, indicating that the localization 
mechanism is intrinsic and arises solely from spin scat- 
tering. In fact, extensive measurements 0, 0, E3 E3 
over the last 10 years indicate that once superconduc- 



tivity is destroyed by the application of a large mag- 
netic field, only two electrically distinct phases exist in 
the cuprates: 1) an insulator with a logarithmically di- 
verging resistivity of the form hiT^/T throughout the 
pseudogap region, x < x c , and 2) a Fermi liquid metal 
for x > x c , where x c w x op t, the optimal doping level. 
In the absence of a field, the most recent [llj in-plane 
transport data on untwinned crystals of YBa2Cu30g.35, 
a composition right at the edge of the superconduct- 
ing dome, corroborate the diverging ln(Xb/T) behaviour 
found in the high magnetic field limit for both p a and p^. 
These authors |llj conclude that the localization mech- 
anism is independent of field and likely to be a conse- 
quence of an intrinsically insulating pseudogap at T = 0. 
Scanning tunneling experiments 12] are consistent with 
the deepening of the pseudogap as the temperature is 
decreased. Finally, recent ARPES measurements^^ 
have detected a finite gap over the entire Brillouin 
zone, including along the d x 2_ y 2 nodal line, in the nor- 
mal state of La2_ K Sr K Cu04, C^-zNazCuC^C^, and 
Nd 2 - a; Ce 2 ;Cu04. In Ca 2 _ a ;Na a ;Cu02Cl2, the gap was 
observed to close at x = 0.12. Similarly, in optimally 
doped Bi2212[3 ) the imaginary part of the self energy 
is momentum dependent but it remains non-zero even 
along the nodal directions. The ubiquity of a complete 
gap in the normal state of both electron and hole-doped 
cuprates prompted Shen, et. al.^j| to conclude that 
gapped excitations in lightly doped Mott insulators is 
a generic feature. 

Quite generally, a pseudog ap liH Hfil | is an example of 
an orthogonality catastrophe 17] . Typically, orthogonal- 
ity leads to vanishing of both the quasiparticle weight, Z, 
as well as the conductivity at T = 0. In this case, we find 
that the orthogonality in a doped Mott insulator arises 
because hole transport is limited by the triplet-singlet 
energy gap. The length scale underlying this energy gap 
is the average separation between doubly occupied sites. 
As this length scale diverges in projected models but re- 
mains finite in the Hubbard model, the pseudogap is tied 
to a non-commutativity of U — > oo and L — ► oo. This 
lack of commutativity offers a possible explanation why 
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all simulations thus far on the t — J model find metallic 
transport 0, 0] near half-filling whereas for the Hub- 
bard model, an insulating state obtains. 

The starting point for our analysis is the Hubbard 
model. In this context, we have been refining |2*H a 
non-perturbative resolvent method^ for calculating the 
single-particle spectral function A(k, uS) = —lmFT(6(t — 
t')({ci a (t), c'j a (t 1 )}) that is based on a self-consistent de- 
termination of the electron self-energy using the Hubbard 
operators. Here, Ci a is the electron annihilation oper- 
ator and FT represents the frequency and momentum 
Fourier transform. In the spirit of cellular methods[23j|. 
the essence of our procedure is to expand the electron 
self-energy for the 2D lattice in terms of the resolvents 
for a small cluster. In our work, the eigenstates of a 
two-site cluster were used to expand the operators in the 
self-energy. As the self-energy can be written as a prod- 
uct of two operators, each of which can be centered on 
different lattice sites, a two-site expansion for each op- 
erator captures local correlations (albeit in a pair-wise 
fashion) over at most four lattice sites. Such a local ex- 
pansion has been shown to yield a heat capacity of the 
ID Hubbard system[24| in excellent agreement with the 
Bethe ansatz as well as a pseudogap|21| in the 2D Hub- 
bard model. Our emphasis here is on using the spectral 
function to calculate the conductivity. To obtain a direct 
link between the conductivity and the spectral function, 
we work with the non-crossing approximation 

Rea xx (0 + i5) = 2ne 2 J d 2 k J du'(2t sin k x ) 2 

to the Kubo formula for the conductivity where f(u>) is 
the Fermi distribution function. Although Eq. © is 
only approximate, as it does not include vertex correc- 
tions, we will show that our conclusions are independent 
of any approximation used to compute the conductiv- 
ity. Shown in Fig. Q is the resultant computation of 
the resistivity as a function of temperature for fillings of 
n = 0.97, 0.95, 0.9, 0.85, and n = 0.8. At high temper- 
atures the resistivity increases algebraically regardless of 
the filling. However, at low temperatues, a divergence in 
accord with Eq. (JTJ for fillings close to n = 1 obtains. 
We will determine the crossover filling on general grounds 
later. In contrast, similar cluster treatments of the spec- 
tral function |is| of the t-J model coupled with Eq. 
find a metallic conductivity at all fillings, even arbitrarily 
close to half-filling. In fact, recent [l9( exact diagonaliza- 
tion calculations on finite clusters confirm the inherent 
metallic behaviour at low temperatures, regardless of fill- 
ing, in the t-J model. Metallic behaviour in the t-J model 
is consistent with the extensive numerical [25L Eil Wl\ and 
self-consistent Born calculations!^ which have found 
that a single hole is mobile in a quantum antiferromagnet 
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FIG. 1: Resistivity as a function of temperature calculated 
according to Eq. @ for the Hubbard model (with U — Wt) 
using the spectral function computed previously by Stanescu 
and Phillips 21] for fillings n = 0.97,0.95,0.9,0.85,0.80. The 
inset indicates that the density of states at the chemical po- 
tential vanishes as the temperature decreases giving rise to an 
insulating state for n ~ 1. po ~ h/e 2 . 



described by the t-J model with a quasi-particle residue 
that scales as Z oc J/t where J = At 2 /U . 

What then is the origin of the insulating state for the 
Hubbard model in the underdoped regime? The inset in 
Fig. ||TJ demonstrates that the density of states at the 
chemical potential plummets to zero exponentially as the 
temperature decreases. The conductivity, Eq. is a 
product of the derivative of the Fermi distribution func- 
tion and the spectral function. Because the former is 
peaked while the latter is zero at the chemical potential, 
the product necessarily vanishes leading to an insulating 
state. This cancellation persists to all orders of perturba- 
tion theory. Hence, the insulating state found here is not 
an artifact of the approximate form of Eq. Pjl: rather it 
arises simply because D(ep) — a,t T — 0. 

Because the electron self-energy is expanded in the 
level operators^ for a two-site cluster, we can de- 
termine which local two-site correlations determine the 
physics of the vanishing of the density of states. The 
solid line in Fig. @ illustrates clearly that the chemical 
potential lies in a local minimum in the single-particle 
density of states. This state of affairs obtains because 
nearest-neighbour singlet states (solid circles) and triplet 
(open squares) contribute to the density of states just 
below and above the chemical potential, respectively as 
shown in Fig. J5J. Because the triplet and singlet are 
split by an energy J = 4t 2 /U, their contributions to the 
density of states cannot occur at the same energy. The 
density of states must have a dip which must constitute 
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FIG. 2: Single-particle density of states when the contribu- 
tion from various the nearest-neighbour singlet, FFa, and 
triplet states, FFs, are isolated. Elimination of the singlet 
contribution (solid circles) diminishes the density of states 
just below the chemical potential, whereas the triplet contri- 
bution (open squares) appears above the chemical potential. 
Since these states have an energy splitting of J = At 2 /U, their 
contributions to the density of states results in a real gap at 
T = O.The occupancy (see inset) in FFa and FFs crosses 
exactly at T* . 



a real gap at T = 0. The inset illustrates that precisely 
at the temperature (see Fig. ©) at which the dip in 
the density of states obtains, the occupancy in the ex- 
cited triplet states drops below that of the singlets. This 
definitively proves that it is the singlet-triplet excitation 
gap that limits hole transport in a doped Mott insulator. 
Such a pseudogap can be thought of as a spin gap0 as 
in the context of a spin liquid|29j. Also consistent with 
our finding here is the ferromagnetic polaron picture [3fj|. 
However, neither experimental nor theoretical|3l| work 
supports the ferromagnetic polaron model in the param- 
eter range of the cuprates. In our simple picture that it 
is the local singlet-triplet splitting that gives rise to the 
pseudogap, we expect the corresponding gap arising from 
the orthogonality to be isotropic in momentum space. As 
illustrated by the inset in Fig. © , the curvature of the 
density of states at the chemical potential is positive at 
each momentum indicating that all momenta contribute 
to the pseudogap, though with differing weights. This is 
consistent with the extensive ARPES study of Shen, et. 
al.^3- I n the context of the cuprates, we propose that 
any anisotropyj^ seen in the pseudogap is absent at T* 
but arises at lower temperatures as a result of any order- 
ing phenomena|33l| or pairing that might supervene on 
the pseudogap phase. In fact, others |3J| have concluded 
recently based on cluster calculations on the Hubbard 
model that a pseudogap arises entirely from local corre- 
lations independently of any ordering or pair formation. 



FIG. 3: Doping dependence of the pseudogap energy scale, 
T*. The dependence obeys the functional form, J(l — 
x^ 2 Q /4a 2 ), where £do is the average distance between dou- 
bly occupied sites, the percentage of which is indicated in 
the table. The inset shows the curvature of the density of 
states at the chemical potential as a function of momentum. 
Positive curvature at all momenta indicates that all momenta 
contribute to the pseudogap. 



Two natural questions that arise from this work are 1) 
why do analogous cluster or exact diagonalization studies 
of the t-J model show no indication of localizationflil fl9| 
and 2) what sets the length scale for the energy gap. 
Both of these questions have the same answer. With- 
out the triplet contribution, the pseudogap in Fig. © 
vanishes. However, the triplet contribution lies above 
the chemical potential and hence ispart of the addition 
spectrum. The addition spectrum[l| of the low-energy 
spectral weight (LESW) is a sum of two distinct pro- 
cesses each involving spectral weight transfer between 
the upper and lower Hubbard bands: 1) a static part 
arising from state counting which grows as 2x but more 
importantly 2) a dynamical part that arises entirely from 
the hybridization. Since the triplet is present only when 
t ^ 0, the triplet contribution to the LESW is purely dy- 
namical. In projected models in which double occupancy 
is eliminated at second order, the LESW scales exactly as 
2x|l|. Hence, the dynamical contribution to the spectral 
weight transfer is absent. However, the dynamical contri- 
bution to the addition part of the LESW can be treated 
perturbatively as first shown by Harris and Laneeplij. 
Perturbation theory alone is insufficient to generate a 
gap in an excitation spectrum since the opening of a gap 
represents a phase transition. The essence of the prob- 
lem is that as long as the insulating state is tied to the 
dynamical contribution to the spectral weight transfer 
between the upper and lower Hubbard bands, the length 
scale, £do, over which transport is governed by double 
occupancy must be finite. That is, the physics is sensi- 
tive to the order of limits of U — > oo and L — > oo. Such 
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non-commutativity signals a breakdown in perturbation 
theory as advocated previously |36j. U — > oo, L — > oo 
results in £do > L, metallic transport. In the reverse or- 
der of limits, £do < L and localization obtains provided 
that the rih£,^ < L 2 , n./, — x(L/a) 2 the number of holes. 
n h£,do ~ L 2 defines the percolation limit. By calculating 
the percentage of doubly occupied sites, we obtained £do 
numerically and plotted the T*-line, J(l — cx(£do/a) 2 ), 
in Fig. |PJ. The agreement of this phcnomcnological fit 
with the resisitivity data in which a metallic state obtains 
at x = 0.1 and the crossing in Fig. lend credence to 
our assertion that £do is the relevant length scale for the 
pseudogap. Finally, the scaling form 



(3) 



for the one-hole quasiparticle weight lays plain that the 
discrepancy between the t— j[2j||23,|28| and Hubbard 
results is one of lack of commutativity. In the t — J model 
(no double occupancy), U — ► oo, L — > oo and Z remains 
finite. In the reverse order of limits (Hubbard model) , Z 
vanishes. 

Indeed, other proposals for hole localization exist. 
Some have ar gued that in the t-J model, a hole creates 
a phase string[37|. However, such an exotic state is not 
borne out by extensive numerical simulations on the t-J 
model [25j|. In the spin-fermion model, selective gapping 
occurs at hot spots indicated by the intersection of the 
Fermi surface arcs with the reduced diamond-shaped AF 
Brillouin zone[3^| whereas in the spin-bag model[3^|, a 
gap occurs only along the (n,n) direction. Neither of 
these models, however, possesses the strong correlations 
intrinsic to the doped Mott state. 

To conclude, our proposal that an orthogonality be- 
tween the singlet and triplet states necessarily requires a 
finite length scale over which transport is governed by 
the distance between double occupancies implies that 
U — ► oo and L — » oo do not commute. The emer- 
gence of such a finite length scale in the transport prop- 
erties offers a possible resolution of the breakdown of the 
one-parameter scaling picture^(| f° r quantum criticality 
in the cuprates. Finally, experiments (4l| demonstrating 
that diamagnetism in the pseudogap phase does not per- 
sist all the way to T* proves that pair fluctuations alone 
cannot account for the pseudogap. As advocated here, 
the pseudogap arises from Mottness and any relationship 
between orderingj^ or pairing and the pseudogap is one 
of supervenience. 
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